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We consider the electric and magnetic field fluctuations in the vacuum state in the region external 
to a half-space filled with a homogeneous non-dissipative dielectric. We discuss an appropriate limit 
to an ideal metal and concentrate our interest on the renormalized field fluctuations, or equivalently 
to renormalized electric and magnetic energy densities, in the proximity of the dielectric-vacuum 
interface. We show that surface divergences of field fluctuations arise at the interface in an appropri- 
ate ideal conductor limit, and that our limiting procedure allows to discuss in detail their structure. 
Field fluctuations close to the surface can be investigated through the retarded Casimir- Polder 
interaction with an appropriate polarizable body. 
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I. INTRODUCTION 



Vacuum fluctuations are an outstanding consequence 
of the quantum theory of the electromagnetic radiation 
field, and observable manifestations of vacuum fluctu- 
ations include radiative level shifts and Casimir forces 
P, Q • These forces are long-range interactions between 
neutral objects due to their interaction with vacuum fluc- 
tuations; they are quantum effects of the radiation field 
and have not a classical analogue. Usually Casimir ef- 
fect refers with this kind of interactions between macro- 
scopic objects [jj, and Casimir-Polder forces with the 
interactions of atoms or molecules with a macroscopic 
object such as a surface or among atoms and molecules 
[H- Vacuum fluctuations are changed by the presence of 
boundary conditions given by dielectric or metallic sur- 
faces. This change of electric and magnetic field fluctu- 
ations gives rise to energy shifts of atoms or molecules, 
placed near the surface, that depend from their position, 
yielding Casimir-Polder forces [5|. Evaluating field fluc- 
tuations around a boundary allows also to obtain how 
interactomic Casimir-Polder forces are modified by the 
presence of the boundary @, Q- These facts give fur- 
ther elements making relevant studying the dependence 
of vacuum fluctuations or vacuum energy densities on the 
distance from a dielectric or conducting surface. 

In the case of a perfectly conducting infinite plate, 
the vacuum electric and magnetic fluctuations (E 2 ) and 
(B 2 ), after subtraction of the homogeneous zero-point 
terms existing even in the absence of the plate, behave 
as z~ 4 , z being the distance from the plate; they do di- 
verge at the metal- vacuum interface. The physical origin 
of this divergence and the possible presence of further sin- 
gular terms at the interface has been recently questioned 
in the literature @, |9| , and possibility of their regulariza- 
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tions has been investigated in the case of a scalar field 
using an appropriate potential to represent the wall [lOj . 
It has been shown that, in the case of a half-space filled 
with a non-dissipative dielectric material characterized 
by a real dielectric constant e independent from the fre- 
quency, such divergences at the interface are still present 
in the limit e — > oo (Tlj . An analogous behavior has 
been obtained also if dispersion is introduced using the 
plasma model [f2j]. In the case of a scalar field, it has 
been shown that the stress-tensor components are regu- 
larized by reflection and transmission coefficients; also, 
the surface divergence associated to a perfectly reflect- 
ing mirror is canceled by a new energy density near the 
surface [Hj|. Surface divergences in the Casimir energies 
has also been considered by introducing extra terms in 
the Lagrangian in the form of a (5-function potential, in 
order to simulate the boundary condition [lj]. Research 
in this direction has been also motivated by the quest 
for situations where negative energy densities occur [ll| , 
and solving discrepancies between total self-energies and 
the local energy densities [1, HH . The presence of surface 
divergences may be relevant also from the point of view 
of the cou plin g of the related field energy densities with 
gravity |,U3. 

In this paper we shall consider electric and magnetic 
field fluctuations in the vacuum region near a half-space 
filled with a homogeneous nondissipative dielectric char- 
acterized by a real refractive index n. The field fluctua- 
tions are evaluated using the Carniglia-Mandel triplets as 
field modes. A time-splitting procedure is used in eval- 
uating the frequency integrals, and we notice that this 
is mathematically equivalent to introduce an exponential 
cut-off function e - ^, with — uj c playing the role of 
a cut-off frequency. We then consider the limit n — > oo, 
leading to a model for a metallic material. If the quantity 
1 / 77 introduced by the time-splitting procedure is kept 
finite, this is representative of a situation where the con- 
tribution of field modes with to > ui c is suppressed by the 
exponential cut-off. We can then consider the limit of an 
ideal metallic plate by taking the limit I/77 = u c — > 00. 
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After this limit, the cut-off function is one for all fre- 
quencies and we recover the nondispersive case. With 
this procedure, the ideal metal case is obtained through 
a limiting procedure, and at no point we need to use 
boundary conditions for ideal metals. We find that in 
the limit n — > oo the renormalizcd field fluctuations are 
finite for any finite value of u> c . We are also able to show 
that in the ideal metal limit lo c — > oo surface divergences 
at the metal- vacuum interface emerge, and to investigate 
their properties. As far as we know, the structure of 
these surface divergences cannot be derived if the field is 
directly quantized with a perfectly conducting boundary. 
Moreover, no discrepancy between total self-energy and 
local energy densities is present in our model. Such dis- 
crepancies indeed occur when energy densities are evalu- 
ated directly for the ideal metallic plate, because of the 
singular behavior at z = 0. From the field energy density 
we can also evaluate the electric and magnetic Casimir- 
Polder force on an atom placed near the dielectric half- 
space. Strong radiative interactions of cesium atoms near 
dielectric boundaries have been recently investigated ex- 
perimentally fl8j . 

In Sec. |TT] we calculate the vacuum fluctuations of 
the electric and magnetic field in the vacuum space in 
presence of a nondissipative homogeneous dielectric half- 
space. We use quantization in terms of the Carniglia- 
Mandel triplets, that include evanescent waves, and we 
use the known time-splitting procedure. In Sec. IIIII we 
consider the ideal metal case by a limit procedure from 

I 



where the index L, R indicates waves propagating toward 
the dielectric- vacuum interface from left or right, respec- 
tively, and j = 1 specifies transverse electric (TE) modes 
and j — 2 transverse magnetic (TM) modes. We have 



the dielectric case, and we analyze in detail the field 
fluctuations (and field energy densities), in particular in 
proximity of the metal- vacuum interface. The emergence 
and structure of surface divergences at the interface in 
the ideal conductor limit is discussed in detail. 



II. QUANTUM FLUCTUATIONS OF THE 
ELECTROMAGNETIC FIELD NEAR THE 
DIELECTRIC- VACUUM INTERFACE 

We consider an half-space filled with an homoge- 
neous dielectric medium and the vacuum in the other 
half-space. We label z the direction orthogonal to the 
dielectric-vacuum interface: z < is the dielectric half- 
space and z > the vacuum half-space. Our model has a 
translational symmetry in the x-y directions. We assume 
that the dielectric is non-dissipative with a (real) dielec- 
tric constant independent on the frequency. We shall 
later discuss how to introduce in our model a dependence 
on the frequency, in particular when considering the limit 
for an ideal conductor. 

First step is the calculation of the square of the electric 
and magnetic fields in the vacuum half-space (z > 0) in 
the ground state of the field. We use the well-known 
quantization scheme in terms of the Carniglia-Mandcl 
triplets [l9j]. The Carniglia-Mandel modes are given by 
the following expressions 
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also defined the polarization operators 
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appropriate polarization unit vectors. 

In order to distinguish waves propagating in the posi- 
tive and negative directions of the z axis, we have defined 
in the free space (z > 0) the wavevectors 



— [k x ; k y ] ±fc z ) — (k|| ; ±fc z ). 
Inside the dielectric (z < 0) the wavevectors are 



with 



kdz =^J(n 2 - l)k 2 + n 2 k 2 z 



~Jk 2 dz -(n*-l)k 2 



(3) 



(4) 



(5) 



At the dielectric- vacuum interface (z = 0), k|| and 



the frequency uik are conserved, while k z appropriately 
changes in order to ensure the required continuity of E||, 
Dj^ and B. For k 2 dz less than (n 2 — l)fc 2 , k z becomes imag- 
inary and a total internal reflection of L-modes occurs at 
the interface yielding an evanescent wave in the vacuum 
region. Thus evanescent waves are correctly taken into 
account by using the Carniglia-Mandel modes. Proof of 
orthogonality and completeness of these modes can be 
found in ■ Similar modes have recently been used 

for the quantization of the electromagnetic field in the 
presence of a nondispersive and nondissipative dielectric 

slab fm, m. 



The explicit expression of the modes for z > (vacuum 
region) are 
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In terms of the modes ([T]) , we can write the expressions of the electric and magnetic field operators 



(6d) 



E(r,t) = V f rf 3 k a kj e"^' fe(kj, r) - oLe 4 "' 1 f£(kj,r) 

j ■ /fc ->° L 

+ ^ / d 3 k d iV27rfej fc L dJ e~^' f L (k d j,r) -ft^.e^*' f£(k d j,r) 
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B(r,i) = J2f d3k 



1 2nhc 2 



a kj e- iWfct V x f fl (kj,r) + aLe to "' t V x f fl( k ^ r ) 



27T?iC 2 



& kdJ e-*-** V x f L (k d j,r) + b\. V x f£(k dJ ,r) 



(8) 



(operators a^j, and &k d j, refer to R- and L-modes, respectively). 
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The field Hamiltonian, after subtraction of the zero-point energy, is given by 

,■ Jk z >0 ' ,• Jk iz >Q 



(9) 



where the annihilation and creation operators for pho- 
tons in the different field modes satisfy the usual bosonic 
commutation rules 



[ b * d j: b l 



3 



5 r y <5 3 (k rf -k^) 



(10) 



(all other commutators vanish). 

Using the expression ([7]) for the electric field with the 
modes (O, we can evaluate the average value of a space- 
time correlation of the electric field in the vacuum space 
near the interface (z > 0) on the ground state of the field, 
obtaining (this calculations is similar to that in and 
we are giving some detail of it in order to show how the 
time-splitting procedure we are going to use is equivalent 
to a high-frequency cut-off) 

(E A (r,*)E A (r',t')) 

= 2nh\Y, I d*k u k e -^(*-0 f> A (kj, r) f>*(kj, r') 
+ W d 3 k d Wfe e -*-*(*-0 f> A (k rf j, r) f> A *(k dJ ,r') 

j Jk dz >0 

(11) 

with A = x, y, z. This quantity diverges for r' — > r and 
t' — » t, but we can use a point-splitting procedure, by 



introducing the following quantity 

(E 2 ) r) = (E A (r, t) E A (r, t' = t + in)) (12) 
with rj > 0. 

In the next section we shall show that this time- 
splitting procedure is mathematically equivalent to in- 
troduce a cut-off frequency in the frequency integrals, 
and this allows us to obtain the ideal conductor case as 
a limit process. Thus, for r\ — > we have (E,),, —5- (E?). 
As long as rj is finite, the integrals do not diverge and we 
obtain 

(E 2 ), = 2** W rf 3 k cj k e-™ |f> A (kj, r) | 2 

j Jk z >0 

+27Th^2 [ d 3 k d u k e-^» |f> A (k dJ ,r)| 2 . 

j Jk d z>0 

(13) 

In this expression, contributions from both traveling and 
evanescent waves are taken into account. 

After some algebraic calculations, using ([5]) and ([6]), 
we obtain the following integral expression for (E 2 )^ for 
z > and arbitrary n 
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An analogous expression can be obtained for the magnetic part (z > and arbitrary n) 
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>L + N 2 *i + » 4 l*»l a _ 

In the next Section we shall consider the appropriate limit of these expressions for the case of an ideal metal. 
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III. THE IDEAL CONDUCTOR LIMIT: SURFACE DIVERGENCES AT THE INTERFACE 



We now consider the limits tt. — 5-1 and n — > oo of (JT4J), respectively yielding the case of vacuum space and of the 
conductor, 



(E 2 )- = lim^E 2 ), = ^ l\ [°dk, fc|| JkfTk* e - c "V^ ; 



JO 



(16) 



ihr r°° r°° h 
(E 2 )™ 11 = lim (E 2 ) r) = — / dh I dk 



o Jo 



k\ + k\ 



(fc 2 + 2fc 2 sin 2 (fc z z)) e~ CI 'V^ 



-k'i 



(17) 



Analogous expressions are obtained from (|15p for the magnetic fluctuations, for the vacuum space and for the 
conductor case, 



(B 2 )™ c 



lim(B 2 )^ = ^ph pk z fc|| ^fcf+fc 2 e - c "v^; 



(18) 



(B 2 )™ 11 = lim <B 2 >, ; = — I dk 



dk z 



k 2 + kl 



(k\ + 2k 2 z cos 2 (k z z)) 



(19) 



The conductor results, both for the electric and magnetic components, can be renormalized by subtracting the 
spatially homogeneous vacuum contribution, that is that obtained respectively in (fT6]) and (|T8| in the absence of the 
material half-space. Thus we obtain 



/xp2\ con 

)r, R 



— [ dku [ dk z 
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(20) 



and 



/■r>2\ con 



/xp2\ con 



(21) 



r 



From the expressions above, it is evident that the 
time-splitting procedure introduces an exponential cut- 
off in the frequency integrals, suppressing contributions 
of modes with frequency u) > lo c = l/rj. If uj c in kept 
finite, the renormalized field fluctuations (|20l) and (|2"Tj) 
are finite for any z, including at the interface z = 0, due 
to the regularization introduced by the cut-off function. 
Also, the electric and magnetic parts are opposite each 
other in all points and for any value of r\. If we take 
r\ = 0, we get a divergence at the conductor- vacuum in- 
terface z — 0: in this cut-off function such as that 
given by the plasma model for a real conductor (assumed 
valid also at very high frequecies), that introduces a uj~ 2 
factor for frequencies larger than the conductor plasma 
frequency, would not be enough to avoid the divergence 
at the interface, as showed in 12]. There is however some 
controversy in the literature about which model of a real 
conductor (plasma, Drude or other) should be used for 
calculating Casimir energies and the Casimir force be- 



tween real materials. 

An explicit evaluation of the integrals in p^|) . (fTT|) . 
pl| yields: 



2\ vac 
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irc 3 r) 4 



for the vacuum case; 



2 \ con 
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I2h ^ 4cft (I2z 2 - cV ) 
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for the conductor case, and 



2 \ con 



(E 2 )- 



4cfi (122 



-V) 



7T (4Z 2 +C 2 r/ 2 ) 3 



(22) 



(23) 



(24) 



for the renormalized conductor case. 

As expected, (E 2 )™ c and (E 2 )™ 11 diverge for rj ->• 0. 

The limit r\ — > of (|24p for z ^ gives back the well- 
known result for the vacuum average value of the electric 
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field squared in the presence of a perfectly conducting 
plate [lH 



(E 2 )?> on = lim(E 2 )™" = — — 



(25) 



Analogous results are obtained for the magnetic field 
fluctuations 



/T>2\vac 12ft 
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'47TZ 4 ' 



(26) 



(27) 



(28) 



(29) 



The renormalized field fluctuations (E 2 )™^ and 
(B 2 )^™ given by ((211) and dHJ), which are proportional to 
the renormalized electric and magnetic energy densities, 
are finite provided 77 is not vanishing. Equations (1251) and 
(l29t show also, as it is already known, that the presence 
of a perfectly conducting plate increases the fluctuations 
of the electric field whereas reduces magnetic field fluctu- 
ations. We can now address the main point of this paper, 
that is the limit to an ideal conductor, and in particu- 
lar the behavior of the field fluctuations near the plate, 
and the presence and origin of surface divergences in our 
limit to the ideal conductor. As already noticed, the use 
of the time-splitting procedure introduces an exponen- 
tial cut-off function in the frequency integrals of SeclTTl 
giving field fluctuations. This is equivalent to assume a 
cut-off function with the scale of the exponential related 
to a frequency such that modes of higher frequency are 
suppressed. A similar cut-off function has already been 
used in [l3| and in [Io[ for the scalar field case. Thus 
we can take expression ((24)) and (l28l) as the renormalized 
fluctuations of the electric and magnetic field for this 
model of a conductor with an exponential cut-off func- 
tion. In this model, I/77 plays a role analogous to that 
of a plasma frequency (but we are not dealing with the 
plasma model, that yields a different form of the cut- 
off function). The fact that an exponential cut-off could 
be not representative of a real conductor is not essential 
in our case, because our interest is to consider the limit 
77 — > or u> c — > go: in this limit our cut-off function 
is approaching one for all frequency, and we recover the 
physical situation of a nondispersive material. However, 
our procedure of considering the ideal conductor through 
a limiting process allows us to show the existence of sur- 
face divergences, and discuss their structure. As we shall 
show in the following, this does not seems possible if the 



field is directly quantized with the perfectly conducting 
plate, nor if the plasma model is used in the limit of a 
large plasma frequency (because, as it is shown in [1^ . 
the plasma model does not completely remove the diver- 
gences at the metal- vacuum interface). 

We first notice that, for any r\ > 0, the integrals of 
the renormalized energy densities over all vacuum space 
vanish for both their electric and magnetic parts, because 



dz 



12z 2 



1 2 
c 77 



(4z 2 + cV) 3 



0. 



(30) 



Thus both electric and magnetic renormalized energy 
densities vanish when integrated over the z > half- 
space, whichever the cut-off frequency is. This also 
means that, in the ideal conductor limit ij — > 0, the (di- 
verging) electric and magnetic energy for z > 0, propor- 
tional to the spatial integrals of (|25p and (T2"9")l respectively 
(both behaving as z -4 for z ~ 0), must be canceled by 
electric and magnetic energies confined at z — in the 
form of surface divergences at the conductor- vacuum in- 
terface. The existence of surface divergences was already 
guessed in [l4| in the case of an ideal flat boundary con- 
dition for a massless scalar field; additional 5— function 
terms were added in the Lagrangian in order to take into 
account the presence of the plate. A similar situation 
occurs for a scalar field with appropriate diverging po- 
tentials which simulate the boundary 10]. Our approach 
of obtaining the ideal conductor through an appropriate 
limit process, as we will now show, actually allows us to 
physically understand the origin and properties of such 
surface divergences of the electric and magnetic energy 
densities, without additional hypothesis. 

We now consider the behavior of the electric fluctua- 
tions and energy density for growing values of I/77, and 
compare it with the same quantity for an ideal metal as 
given by ((23)) . 

Figured] shows (E 2 )Jj™ as given by (|2~4")) with a value of 
1 jr\ comparable with a typical plasma frequency of a real 
metal, in comparison with the z~ 4 behavior of (E 2 )^ 311 
given by (j25[) . Significant differences are evident in the 
proximity of the interface. 

Figure[5]shows the behavior of (E 2 )^™ for different val- 
ues of 77, and compared with (E 2 )p on (continuous line). 
For any nonvanishing value of 77, i.e. for any finite value of 
the cut-ff frequency w c , renormalized fluctuations and en- 
ergy densities are finite in all points of the vacuum region 
and there are not divergences at the surface 



a maximum at z™ 8 ^ 



IjC 

~2 

v min 



e2 )™r. has 

ch 



~Tr~ with the positive value — -rr 



and a minimum at z™ ln = with the negative value 
— — -tt. The width of the curve, that can be estimated 
as the distance between the two inflection points around 
the maximum, is A, ; ~ 0.5?7C. For an increasing cut- 
off frequency I/77, the curves in figure approach the ideal 
conductor limit for large distances from the interface, but 
significant differences still remain close to the surface. In 
fact, while the ideal conductor limit (E 2 )p on diverges with 
positive values at the surface, (E 2 )™^ assumes more and 
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FIG. 1: (Color online) Comparison of <E 2 )™r with I/77 = 
2xl0 16 Hz (dashed blue line) and (E 2 )^ 11 (continuous black 
line). 



/ 172 \ con 
),,R 




l/i]= 1.2xl0 16 Hz 
I/j/= 1.6xl0 16 Hz 
1/;j = 2.0x10 16 Hz 
l/^ = 2.4xl0 16 Hz 
l/i/ = 2.8X10 16 Hz 
l/i/ = 3.2X10 16 Hz 
1/?; -> 00 

z [nm] 



FIG. 2: (Color online) Comparison of (E 2 )=°£, for different 
values of 77. When 77 decreases (cj c increases), renormalized 
field fluctuations and energy densities for z > tend to the 
z~ 
at 



4 law of the ideal conductor case (continuous line), whereas 
~ surface divergences appear with a positive and a 



negative peak that in the limit tend to squeeze at 



0. 



more negative values as 1/r] — > 00, and the width of the 
curve reduces to zero. Maximum and minimum values 
of (E 2 )™^ tend to collapse each other at the surface in 
the ideal conductor limit, yielding a surface divergence 
containing a nonvanishing electric and magnetic energy 



-> z" 



0. 



(31) 



Thus, by analyzing the limiting case ui c — > 00 (77 — > 0) it 
becomes evident that the well-known diverging behavior 
as z~ 4 near the surface for the ideal conductor, indeed 
originates from the maximum of (E 2 )™^. This is clearly 
shown in Figure [5] The negative divergence at z = 
when 77 — » is completely lost if the calculation is di- 
rectly performed for an ideal conductor. Our approach 
thus makes clear the origin of the surface divergences of 
the renormalized squares of electric and magnetic field; 
these divergences, and their detailed structure, naturally 



appear in the limit process from the dielectric to the ideal 
conductor that we have used. 

The limit lj c — > 00 is well defined in our model for any 
distance z except z = 0. Also, due to (j3T))) . both integrals 
over all vacuum space of the (renormalized) electric and 
magnetic energy densities vanish, as expected from the 
evaluation of global field energies. This feature does not 
seem to occur if the calculation is performed directly for 
an ideal boundary or the plasma model is used even in the 
limit of an infinite plasma frequency (see [HI), because in 
these cases the surface divergences are not fully included 
in the expressions of the field energy densities. By taking 
into account the surface divergences with our limit pro- 
cedure, we thus obtain consistency between global field 
energies and (integrated) local energy densities. 

Field fluctuations and energy densities near the surface 
can be investigated through the retarded Casimir-Poldcr 
interaction energy with an appropriate polarizable body 
with static polarizability a. When the distance d between 
the surface and this body is larger that the wavelength 
associated to its main transition, the so-called far zone 
approximation holds; in this case the Casimir-Polder en- 
ergy is given by AEe = — a-(E 2 )^ 01 r \ (d)/2 (see for exam- 
ple [HI). A similar expression holds for the magnetic 
fluctuations. 



IV. CONCLUSIONS 

In this paper we have considered zero-point electric 
and magnetic field fluctuations (or equivalently field en- 
ergy densities) in the vacuum space, when half-space is 
filled with a homogeneous non-dissipative dielectric. The 
material is characterized by a constant real refractive in- 
dex n. We have then taken the limit n — > 00 and, by 
introducing a time-splitting procedure with a parameter 
77, we have mathematically included also a high-frequency 
exponential cut-off function characterized by a cut-off fre- 
quency U3 C = 1 jr\. We have considered in detail two suc- 
cessive limits: n — ► 00 with ui c finite, and then uj c — » 00 
(ideal conductor limit). We have found that no diver- 
gences in the renormalized field fluctuations (and in the 
renormalized field energy densities) are present if uj c is 
kept finite. In this case, field fluctuations at small dis- 
tances from the interface significantly differ compared to 
known results for a perfectly conducting plate, as shown 
in Figures [T] and [2] In the limit cu c — > 00 (ideal conductor 
limit) surface divergences of the field fluctuations (or en- 
ergy densities) at the interface z = are found, while for 
z > the fluctuations approach the well-known be- 
havior of the perfect conductor case. The structure of the 
surface divergences has been discussed in detail. These 
surface divergences of field fluctuations are not obtained 
if the field is directly quantized in the presence of an 
ideal metallic surface; also the plasma model does not 
seem to allow to obtain their structure. Thus our ap- 
proach of obtaining the ideal metal through an appropri- 
ate limit procedure starting from a dielectric has allowed 
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us to obtain in a natural way the explicit structure of the 
surface divergence at the interface between vacuum and 
a conducting material. Having taken into account the 
field energy in the surface divergences has also allowed 
us to show consistence between global and local field en- 
ergies. Finally, we have stressed that field fluctuations 
and energy densities near the interface can in principle 
be investigated through the Casimir-Polder interaction 
energy with an appropriate electrically or magnetically 
polarizable body placed near the interface. 
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